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THE  STABIUn  OF  OPERATION  OF 
A  LIQUID  FUEL  ROCKET  ENQip 


/^he  foi^lowing  le  a  translation  of  •Uat^hlvoat*  rabot/ 
ahidkoatnogo  raketnpgo  dTlgatoIya*  (Etaglish  vara Ion 
above)  hy  K;I.  Artononov  in  la  AN  SSSR.  QtdalflBlTO 

XsJtliaJSBuki  tigManlto  1  ^  if 

Feb,  1961,  pagep.  64-69.^ 

The  liquid  roclcet  engine,  (LRS)  is  eoneidered  as  an  osoill- 
atory  systen  with  delayed  feedback;  toahsient  processes  in  one 
of  the  eleamits  of  the  engine  (the  fuel  lines)  are  described  by 
neans  of  a  wave  equation.  The  lag  and  the  wave  processes  induce  an 
alternation  of  stable  and  unstable  regies,  as  the  thrust  of  the  UlE 
is  reduced.  The  oscillation  frequencies  are  then  close  to  the  natur¬ 
al  i^equencies  of  the  fuel  pipelins.  If  the  lag'  is  a  li^ge  ope,  it 
is  only  possible  to  have  relatively  low-frequency  oscillatory  rsgiines 


The  self -excitation  of  oscillations  in  an  IRE  is  to  a  large 
extent  detenoined  by  the  diaracter  of  the  process  of  converting 
liquid  fuel  .c6]q>onents  to  gaseous  combustion  products.  Before 
being  COTvwted  to  gaseous  (Hmbustion  products  the  oxidlser  and 
the  fuel  must  go  through  a  series  of  preparatory  processes  (aton- 
iaation  of  the  liquid  jets,  heating  and  evaporation  of  the  droplets, 
mixing  of  the  components,  chemical  reactions,  etc,).  The  volume  of 
fuel,  passing  through  this  stage,  has  practically  no  effect  on  the 


g«ft«ou«  prodttotif  tlM  Inflow  of  whieh  olao  dotoxvinaa  iha  ohaabar 
praaauxw*  Glaan  suffleiantly  rapid  Ttrlationo  la  fuel  eofWvaiptioD 
this  naans  that  at  any  glaaa  nonant  tha  ohaabar  praaaura  la  datana> 
laad  by  tha  eoaatiaptioa  of  fual  during  a  oartaln  praeadlng  Intanral 
of  tlaa*  In  othar  words,  tha  eonnaotlon  faatwaaa  tha  flow  of  fual 
through  tha  atoalsar  and  tha  diaabar  praaaura  la  eharaetarlaad  by  a 
daisy*  In  axlatlng  l£t  fual  food  aystaaa  this  eonnaotlon  is  also 
two-irag^r  fluotuatleos  In  tha  ohaabar  praaaura  also  affaot  tha  rata 
of  flow  of  fual  Into  tha  ohaabar*  Thus,  with  raapaot  to  tha  nature 
of  tha  working  prooaaa  In  tha  ocaboation  ohaabar  and  tha  oharaetar- 
latiea  of  tha  fual  faad  syatan  tha  IRS  la  a  ayatan  with  dalayad 
faadbaek. 

Qna  of  tha  alMoants  of  tha  propulsion  plant  -  tha  fual  Unaa  - 
oc^it  to  ba  oonalderad  an  alaaMint  with  dlstrlbutad  parawatara,  alnea 
tha  transit  tins  for  the  pressure  ware  throu^  tha  plpallne  and  the 
period  of  tha  osoilletlona  are  ooMM&aurable. 

Tha  lag  and  the  vase  proeaaaaa  datandna  a  tarlaa  of  prop- 
ertlaa  of  tha  dynaades  of  an  Utt  whieh  wa  shall  diaeusa  balew  with 
rafaranoa  to  the  sli^laat  type  of  propulaion  plant* 


In  ordar  to  darlwa  tha  stability  oondltlona*  wa  shall  writs 
out  tha  aquations  daaorlbing  non-ataady  prooaaaea  In  an  UUE*  la  is 
U8\;ial,  in  laTeatlgatlng  stability  In  tha  aaall  wa  shall  use  llnaar- 
laad  aquations* 


Wa  shall  eharaetariaa  tha  prooaas  of  eonuartlng  liquid  fuel  to  gas 
fay  eaans  of  tha  ocnearslon  curra  ^(t)  -  ratio  of  tha  aaaa  of  tha 
gas  fomad  at  a  nowent  t  fron  the  ccaibustion  of  a  giran  woluae  of 
fuel  to  the  Initial  naas  of  tha  fuel  introduced  at  tha  aonant  t  •  0. 

£rldantly,  tha  eonvarslon  ounra  ollnbs  anoethly  fron  saro  and 
tends  asynptotloally  to  unity  CPlg*  1)*  This  ounra  la  only  an 
araraged  cbaraotaristle  of  tha  prooessaa  at  work}  Its  fom  la 
datamlnad  both  by  tha  law  of  eonwaraion  of  tha  indiwldiMd  psrtiele 
(tha  drop)  of  fuel  and  by  diffaranoaa  In  tha  oondltiona  aasoclatad 
with  the  eonrarsioQ  of  the  different  partiolaS  (a*g.  due  to  the 
non-uniform  also  of  the  dropa)*  Wa  shall  assuna  that  tha  ooniTeraloo 
curve  is  datamlnad  only  by  ataady-stato  parametara  and  la  not 
affected  by  small  fluctuations* 

If  wa  Introduce  tha  rata  of  ooofveralon  ^  «  db/dl  $  than 
tha  mass  of  gaa  eonvartad  par  unit  of  tiate  il#r(/)  nay  ba  eacpraaaad 
as: 

00 

^r(0-  5  A/♦(^  (2.1) 


vh«r«  M^Ct)  l8  iiiB  uBis  rata  of  InfloM  of  fuel  through  tha 


•pray 


ha«^. 


U  t. 
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Ua  ahall  latroduoa  tha  aaall  dlnanalmilaaa  darlationa 


nff  If  ».• 
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Tha  tndax  0  danotaa  stationary  valuas  of  tha  functions • 
Taldag  into  aoooont  that 


to 


va.  gat  Instaad  of  (2.1) 


CO 


(2.2) 


We  shall  eonsldar  that  tha  gas  porassura  is  uniform  throughout  tha 
ehaabar  at  all  times.  This  is  oorraet  for  osolllationsy  tha  period 
of  which  is  muoh  greater  than  the  time  tha  pressure  wave  takas  to 
traverse  the  length  of  tha  chamber. 

The  equaticm  of  conservati^i  of  mass  of  the  gas  in  the  chambe: 
has  the  form 


+  j»/e  =  il/r 


(2.3) 


where  Qg  is  the  mass  of  gas  in  the  chamber  and  Ks  is  the  mass  rate 
of  flow  of  gases  through  the  aupersonio  nosale. 

If  during  the  oscillation  the  relationship  between  the  fuel 
ooaqponents  does  not  change,  it  is  possible  to  consider  the  ohamber 


prooMsos  lsoentropic«  i*6« 


»  const 


(2.4) 


To  equations  (2«4)  and  (2.5)  wa  auat  add  the  equation  of  atataj 


and  Mpresaiotts  for  and  Kq 


r~l 


(2.5) 


where  pK>  fx  and  are  reapeotiweljr  the  pressure,  density  and 
teoperature  of  the  gases  in  the  e<»bustloa  dhaabar;  R  is  the  gas 
constant,  ^  the  adiabatic  coefficient,  the  woIum  of  the  coa- 
bustion  duusber,  Fgp  the  area  of  the  critical  cross  section  of  the 
noasle  and  g  the  aeoaleration  due  to  gravity. 

BquatiMis  (2«3)-(2«5)  give  us  the  oonneetlon  between  Um  and 
the  dlnensionlees  pressure  varlat^t  in  the  chaabar 

in  th.  following  fomt 


(2.6) 


Hare  tg  is  the  tlas  the  gases  reaain  in  the  combustion  ohaabat 
divided  tgr  the  adiabatic  coaffieiaiit.  ) 


Sauation  of  fuel  feed  svatem 

It  is  assuBwd  that  the  flew  of  liquid  is  not  influenced  by  inartla 
foroes  due  to  the  acceleration  the  rocket.  We  shall  consider 
a  coapressed  gas  fuel  feed  system  consisting  of  a  fuel  tank,  a 
boaogeneous  pipeline  and  a  noa.-alastie  spray  head. 

The  action  of  the  liquid  in  the  straight  cylindrical  pipe 
is  deteraiaed  by  Zhuleov'a  equations: 


(2.8) 


^  these  equations  f  is  the  oonstant  density  of  the  liquid, 
s  is  the  speed  of  sound  in  the  liquid  in  an  elastic  pipe. 

The  end  of  the  pipeline  at  the  tank  nay  be  considered  closed; 
during  oscillstions  the  preasure  in  this  saotico  is  cctutant: 

X « 0,  *£«>()  (-‘^) 


Neglecting  the  inertle  of  the  liquid  in  the  sprey  head  cavity, 
we  get  a  second  boundary  condition  from  the  expression  for  the  flow 
of  fuel  through  the  spray  head; 


Introducing  the  variables 


(2.  MM 


p«% 


Pm 


and  eliminating  from  (2*9),  we  get  for  m  the  wave  equation 


iPm  -2  3^  j, 
with  the  boundary  conditions t 


(2.11) 


(2.12) 


where 


The  second  condition  in  (2*12)  is  obtained  by  linearising 
equation  (2«10). 


We  shall  apply  the  Laplace  transformation  to  the  equations 
thus  obtained!  retaining  the  previous  notation  for  the  transformed 
functions*  The  operator  for  differentiation  with  respect  to  time 
will  be  designated  q. 

The  wave  equation  (2*11}  then  goes  over  into  an  ordinary  dlff 
erentlal  equation  ,  the  solution  of  which  for  rnCq,!),  equal  to 
taking  into  account  boundary  conditions  (2»12),  is  obtained  in  the 
form:  ^ 

"•* - ^ 

*♦+*.*»?-  ■* 


Instead  of  (2,1)  and  (2,7)  we  get: 

mr  fp  (^) 

m- 

'^-iETK 


(3.2) 

(3.3) 


From  eqa,  (3«1)«‘(3*3)  it  is  easy  to  obtain  the  charaoteristio 
equation  of  the  (qrstem: 


(3.4) 
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4.  Ll«l^  af  at>hl»  ooTBtloB  of  aa  IBM 

Wtt  obtain  the  aquatlona  for  eetabllshlng  the  llalta  of  the 
regions  of  stable  operation  of  an  lAB  and  for  computing  the 
oaoillation  frequency  au  at  the  linlta  of  stability  ty  eubatitutlng 
in  (3.4)  q  s=ri«  - 

These  equations  oan  be  reduced  to  the  form 


Here 


V 


Ijr  0  =  — 


(a 


(6»)  rr  (/a»)  4-  Im»  ^  (iVe).  led- 


(4.1) 

(4.2)  I 

\ 


The  frequency  functions  S(a»)  and  0(w)  are  the  aaplitude 
and  phase  spectra  of  the  curve  ^(t)  respectively. 

To  continue  the  analysis  it  is  necessary  to  assign  a  concrete 
form  to  the  conversion  cunra.  In  many  instances  the  5^(t)  curve 
can  be  epproximatad  by  means  of  a  step  function  (curve  1  in  Fig.  1): 

(p«l,  l>/„ 


This  means  that  the  fuel  la  converted  Instantaneously  into 
gas  tn  seconds  after  admission  into  the  oonbuatlon  ohairdoer.  For 
such  a  curve: 

6’(<0)-fI,  0=-CDt„ 

Approximating  ^(t)  with  function  2  (Fig.  1)  we  get 


For  curve  2,  aa  for  any  aymmetrieal  l|^(t)  curve, 
where  tn  is  the  abaclsaa  of  the  axis  of  ayonstry. 

Let  us  consider  in  greater  detail  the  case  of  a  step  oonvers- 
itm  curve.  In  this  case,  on  Introducing  the  quantity  Z/a  aa  a 
time  scale,  we  shall  have  instead  of  (4«1}  and  (4.2}: 


I.  — r - - 

vtg*a 


(4.3) 


.  _ 

sign  sin  at,,  «  sign  (/^xTk  —  XA,  tg  a) 

lii  ^  *n^  tJ  N 

2r  •  ^n«— »  r^^~} 


(4.4) 


An  folIowB  fron  (4*3  )  and  (4«4)»  tha  atebiliiy  of  the  oper¬ 
ating  regljM  of  an  UUS  is  detenBlned  ^  the  following  four  dimension¬ 
less  parameters: 

<„/  2Ap*  (MV* 

Aa  =  -^,  ^.*7- 

«  <*  ^KO  "ho 

Where  Tfi  ts  the  ratio  of  the  time  lag  to  the  tine  of  peaaege 
of  a  sotind  wave  through  the  length  of  the  pipeline  2 /a; 
ratio  of  the  time  the  gases  stay  in  the  ohiotfber  traduced  'f  times) 
to  2 /a;  h  f  is  double  the  ratio  of  the  pressure  drop  at  the 
spraj  head  to  the  chamber  pressure  and  hg  is  the  i]q«et  pressure^ 
due  to  Instantaneous  total  retaz^ti(»  of  the  liquid,  to  the  chamber 
pressure.  /:  =  /:« //:«"»  pSoS/ 2 would  be  a  more  graphic  hjdraullc 
parameter,  but  in  the  eaalirsis  of  atebillty  it  is  more  ooatrenient  to 
use  the  parameter  hg  . 

Another  parmseter  is  ttie  qtaantity  X  «»«  (7  +  l)/27,  however, 
it  Is  sttff lolently  oorreot  to  put  it  equal  to  unity. 

In  oontrolllx^  ths  thrust  of  an  lAE  by  rsdueing  tbs  pressure 
In  the  fuel  tanks,  the  parameters  Ta  and  hp  remain  e<niBtant, 
the  parameter  h^  is  recced  in  proportion  to  the  pressure  in  the 
ohember  (and  the  thrust);  the  variation  may  be  very  direrse 

and  depends  on  the  type  of  fuel  and  the  foarm  of  atcmlsetlon.  It  is 
possible  to  reckon  that  for  the  msjorlty  of  engines  ths  oontrol  time 
lag  increases  with  the  decrease  in  thrust.  Thus,  it  is  eaqpedient  to 
plot  the  limits  of  stability  in  the  plane  A/ 

Since  hm  is  a  real  quantity,  aquation  (4.3)  gives  the 
limitatlma  on  ^  oscillatlan  frequen^:  the  posslbls  frequenoies 
are  either  close  to  ai«o  <«*  oloae  to  m  »»!!(»«*•  1,  3,...).  T‘or 
these  frequenoies  the  reactive  resistance  of  the  oolumn  of  liquid 
(the  quantity  hg  tan  M  )  is  elose  to  sero.  It  Is  clear  that 

ot »  nn  oorreapoada  to  the  natural  frequencies  of  the  ecoustle- 
al  oscillations  of  ^e  liquid  in  the  pipeline,  i.e.  to  the  frequeneie 

/  =  «35*  (o/s) 

The  value  n  s  1  corresponds  to  the  fundamental,  n  s  2,3,. •• 
to  the  second,  third,  eto.  barm(»io8.  For  low  frsquenoies  the 
function 


A,lga«tiii  (/i 

•nd  thA  ioArtla  of  tho  liquid  «r«  d«fla«d  Iqr  •  slnglo  tiuB  ooMtiult. 
With  rjffleloat  uoourtogr 

lliai&iai/a  fttt  a^Vf^ 

accordingly,  the  oMproaalblllty  of  tbo  liquid  can  be  neglaotod  }xp  to 

3 

By  doteralning  b4  for  a  given  value  of  K  froa  aq*  (4<3)* 
we  oan  find  the  value  of  at  the  llait  of  etabllity  froa  (4*4}  •  i 
The  quantity  la  detendned  correct  to  21-31  (4  1.  l*e« 

It  la  posaible  to  have  oaclllatlona  with  the  sane  dlnenalonleaa 
frequent  ®<  for  several  values  of  ,  differing  from  one  another 
hy  2n/ct.  , 

Fig.  2  shows  the  eonplete  Unit  of  stablUty  In  the  plane  of 
the  paraaetera  Ta  for  h%^»2  and  ta 1  Tor  the 

values  n  »  1,2  and  k  •  1,2.  ^  stability  lies  to  the 

rl^t  and  is  bounded  by  a  aaooth  i^per  curve  (low  frequencies)  and  a 
saw-toothed  curve,  for  which  the  oscillation  frequencies  are  close 
to  the  natural  frequencies  of  the  pipeline.  In  this  particular  case 
the  "teeth*,  corresponding  to  the  valnea  n  *  1,  k  >1  and  n  s  2, 
k  >  2,  do  not  change  the  region  of  stability. 


In  the  o«0e  of  this  psrtioular  aagine,  wh«n  the  preesure  in 
the  fuel  tanks  le  reduced,  the  working  point  (^,  t„)  1»  displaced 
tqiwards  and  to  the  left.  Accordingly,  as  the  t&uat  of  the  IRE  Is 
reduced,  stable  and  unstable  reglaes  nay  alternate ^  An  increase  In 
will  not  take  the  engine  out  of  the  region  of  etabillty,  if  the 
working  point  crosaee  the  i^per  saooth  boundary,  tdiieh  retreats  to 
infinity  as  h^f  approachea  1  (aore  accurately  X  )• 

The  effect  of  other  paraaetcsra  on  the  region  of  stablli^ 
oan  be  eonreniently  traced  by  considering  separately  the  upper  bound” 
ary  and  the  tooth  with  k  «  0,  n  »  1.  The  effect  of  and  hg 
My  be  seen  in  Figs.  2  and  3«  An  increase  in  hg  leads  to  a  narrouin 
of  the  teeth,  without  affecting  their  position  end  length,  and  sharp 
ly  displaces  the  upper  boundary,  at  the  sane  tine  Increasing  the 
nifflber  of  teeth  (value  of  k).  The  jiffect  of  J/a  is  illustrated  in 
the  plane  h^  ,tn  in  Fig.  4«.  The  effect  of  I  /a  on  the  operating 
stability  of  an  IRE  nay  be  very  divarae.  An  increase  in  i  /k  leads 
to  a  reduction  In  the  region  of  stability  due  to  elongation  of  the 
j^th  ai^  powerfully  broadens  it  by  dispacing  the  upper  boundary, 
when  //a  is.  increased,  the  points  of  thie  teeth  ascend,  but  at  . 
large  values  of  ^ /a  it  is  necessary  to  take  into  account  teeth 
corresponding  to  hamonics  of  the  fundaMntal.  It  is  Interesting  to 
note  that,  as  / /a  varies,  f or  n  »  1  the  points  of  tlie  teeth  are  dis¬ 
placed  aloM  a  curve  corresponding  to  2  «  0  (broken  line  in  Pig.  4) 
This  ease  (2*0,  step  conversion  curve)  has  already  been  dis¬ 
cussed  by  M.S.Patanstti. 

As  the  preceding  discussion  shows,  it  is  possible  to  represent 
an  IRE  by  neans  of  the  following  block  diagram: 


liii-S’H-il 

i _ :!iJ 


where  the  links  M,  P  and  X  denote  respectively  the  fuel  suj^ly 
system,  the  conversion  process  and  the  combustion  (dumber  as  a  gas 
capacity.  The  characteristic  equation  of  such  a  ^rstem  can  be 
written  out  at  once  from  the  known  transfer  functions  of  the  links: 

A’  {q)  M  {q)  n  (7)  1 

For  the  propulsion  unit  discussed  above  equations  (3.1), 
(3.2)  snd  (3.3)  give  , 

w (</)=. <r("/). 
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